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ABSTRACT 

ABSTRACT: The COBE detection of microwave anisotropics provides the best way of fixing the amplitude 
of cosmological fiuctuations on the largest scales. This normalization is usually given for an re = 1 spectrum, 
including only the anisotropy caused by the Sachs- Wolfe eff'ect. This is certainly not a good approximation for 
a model containing any reasonable amount of baryonic matter. In fact, even tilted Sachs- Wolfe spectra are not 
a good fit to models like Cold Dark Matter (CDM). Here we normalize standard CDM (sCDM) to the 2-year 
COBE data, and quote the best amplitude in terms of the conventionally used measures of power. We also give 
normalizations for some specific variants of this standard model, and we indicate how the normalization depends 
on the assumed values of re, Qb and Hq- For sCDM we find {Q) = 19.9±1.5/iK, corresponding to Ug = 1.34 ±0.10, 
with the normalization at large scales being B = (8.16 ± 1.04) x 10®(/i~^Mpc)'*, and other numbers given in the 
Table. The measured rms temperature fiuctuation smoothed on 10° is a little low relative to this normalization. 
This is mainly due to the low quadrupole in the data: when the quadrupole is removed, the measured value of 
it(10°) is quite consistent with the best-fitting {Q). The use of {Q) should be preferred over it(10°), when its 
value can be determined for a particular theory, since it makes full use of the data. 
Subject headings: cosmic microwave background — large-scale structure of universe 



^ INTRODUCTION 

During the 1980's, it was standard practice to normalize cos- 
^^ological models at a scale of ~ 10/i~^Mpc, using a quantity 
llated to the clustering of galaxies (here the Hubble con- 
lant Hq = 100/ikms~^ Mpc~^). Due to processing of the 
rimordial spectrum, this method of normalization requires 
S^sumptions about both the equation of state for matter in- 
iJsJde the horizon and the relationship between the observed 
©Hucture and the underlying mass in the universe. After 
Jtke COBE DMR detection of Cosmic Microwave Background 
}2pMB) anisotropics (Smoot et al. 1992), it became possible to 
"^rectly normalize the potential fiuctuations at near-horizon 
^ales, circumventing both of these problems with the 'conven- 
tional' normalization. To be specific, diff'erent authors have 
chosen to use one or other of the following COBE-derived quan- 
tities: (1) it(10°), the rms temperature fiuctuation averaged 
over a 10° FWHM beam; (2) (3rms-ps(ra = 1), the best-fitting 
amplitude for an re = 1 Harrison-Zel'dovich (HZ) spectrum, 
quoted at the quadrupole. We will henceforth refer to this 
quadrupole expectation value as '(<?)'. 

There are reasons for preferring one of these quantities over 
the other. First, cr(10°) has the advantage of not depending 
on the assumption of a model in the fitting of the data, since 
it is purely an observationally-determined quantity. However, 
because of this, it does not provide the most accurate normal- 
ization for any specific theory. Furthermore, one must take 
great care to account properly for cosmic variance before us- 
ing cr(10°) to normalize a model (White et al. 1993). On the 
other hand, {Q) (re=l) gives the best amplitude for a pure HZ 
spectrum of CMB fiuctuations, including only the Sachs- Wolfe 
(SW) effect (Sachs & Wolfe 1967). Although the SW effect is 
the main contributor to large-scale CMB anisotropy, other ef- 
fects are not generally negligible, even on COBE scales. In par- 
ticular, as shown in Figure 1, the sCDM model, with Qb = 0.05 
is not well approximated by a S W spectrum for any re. In other 
words, the low £ multipoles are not purely the SW potential 
fiuctuations even on the angular scales probed by COBE. 
Analysis of the first year of COBE data gave it(10°) = 
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30±5/iK and (Q) (re=l) = 17±5/iK (Smoot et al. 1992, Wright 
et al. 1994a); after the second year of data these became 
(7(10°) = 30.5 ±2.7/iK and (Q)(re=l) = 19.9±1.6/iK (Bennett 
et al. 1994, Gorski et al. 1994a, Wright et al. 1994b). With 
twice as much data there has been significant improvement in 
the error bars. But, confusingly, the best value for it(10°) has 
remained essentially unchanged, while {Q) has increased sig- 
nificantly. A simple calculation shows that it(10°)/ {Q) ought 
to be about 2 for re = 1, whereas the actual ratio is about 1.5. 
There is therefore a discrepancy of as much as 30% between 
these two choices of normalization, as also discussed by Banday 
et al. (1994). 

This apparent inconsistency has led to some confusion about 
the exact normalization to use for a given theoretical model. 
This has implications for the question of how specific models fit 
the large-scale structure data (e.g. Efstathiou et al. 1992), what 
the required bias of galaxies is relative to dark matter, and 
whether particular models are consistent with smaller-scale 
CMB data (e.g. Gorski et al. 1993, Dodelson & Jubas 1993, 
Srednicki et al. 1993, Bunn et al. 1994). It therefore seems 
important to normalize with the full accuracy available from 
the COBE data. 

Other authors have tended to use the COBE data to fit Q 
as a function of re using the SW formula (e.g. Scaramella & 
Vittorio 1993, Seljak & Bertschinger 1993). One approxima- 
tion for standard CDM is to take some value of {Q) for re > 1, 
say re = 1.15 (Bond 1993). However, even this is not an ade- 
quate representation of the sCDM model at small multipoles 
(Fig. 1). The best way to normalize a given model is to carry 
out a detailed fit to the data. The point of this letter is to 
perform this task for the totally 'vanilla-fiavored' CDM, and 
for a range of minor deviations from it, so that it will be clear 
exactly what normalization to use in the future. 

Cosmological Models 

The cold dark matter model (see e.g. Ostriker 1993, Liddle 
& Lyth 1993) has become the 'straw man' model of structure 
formation. In this model f2o = 1; with a variable fraction Qb 
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residing in baryons and the rest in massive (non-relativistic) 
dark matter. The initial fluctuations are assumed to be Gaus- 
sian distributed, adiabatic, scalar density fluctuations with 
an HZ spectrum on large scales, i.e. P{k) oc k"' with n = 1. 
The 'standard' CDM model has come to mean the one with 
Hq = 50kms~^Mpc~^ , as a compromise between age and Hq 
constraints, and with Qb = 0.05, for consistency with the best 
central value from big-bang nucleosynthesis studies (Walker et 
al. 1991, Smith et al. 1993). 
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Fig. 1. — The solid lines are the CMB anisotropy multipoles for CDM 
models for Qb = 0.01, 0.05 ('sCDM') and 0.10 (from bottom to top). The 
dashed lines are SW spectra for n — 1.00 to n = 1.25 in steps of 0.05, as 
computed using equation (4). The curves have all been normalized at Cg, 
which is approximately the 'pivot' point for the COBE data (Gorski et 
al. 1994a). Note that the CDM spectra are not well-represented by any n. 

Although it seems generally successful, one does not expect 
the sCDM model to be precisely the correct description of the 
universe, so it is worthwhile to consider a range of deviations 
from it. First one can allow for variations in the cosmological 
parameters h and Qb- On the scales probed by COBE the 
power spectrum is relatively insensitive to these changes (see 
below). Furthermore, the sCDM model assumes a 'flat' or 
HZ [n = 1) spectrum of density fluctuations. Since inflation 
generically predicts departures from the simple re = 1 form, we 
will quote normalizations for a range of re. 

Should primordial gravitational waves (tensors) turn out to 
be important as well as the scalar fluctuations, their effect will 
be to lower the inferred by roughly the fraction they con- 
tribute to C2, i.e. Cj/C^f. This fraction is currently totally 
unknown, although it may be related to re for speciflc mod- 
els. For example, in 'extended' inflation it is roughly 7(1 — re) 
(Crittenden et al. 1993), although in 'natural' inflation it is 
negligible even if re < 1 (Adams et al. 1992). Since including a 
tensor component also changes the shape of the C^'s, we per- 
form a flt to the data. The two speciflc possibilities mentioned 
above span the interesting range and we quote the normaliza- 
tion in both, cases. 

The addition of a component of Hot Dark Matter will cause 
a negligible effect at COBE scales, but will change the relevant 
normalization at cluster and galaxy scales. Models that have 
a cosmological constant and/or are open (including Baryonic 
Dark Matter models), as well as models containing topologi- 
cal defects, present additional difficulties, and we do not treat 
them in this letter. 

Fitting CDM to the COBE data 

As usual in CMB studies we expand the temperature fluc- 
tuations AT /T[d,(fi) in spherical harmonics J2im'^tmY(,mi9,4')j 
and work in terms of the multipole moments a^m- In general 



the predictions of a theory are expressed in terms of predic- 
tions for the ag^rn- Using rotational symmetry {al^a^im') = 
C(,6in6m'mt where the angled brackets represent an ensemble 
average. If the fluctuations are Gaussian, the predictions are 
fully specifled by giving these C^'s. We have calculated the C^'s 
using a Boltzmann code, and have compared them to other cal- 
culations (N. Sugiyama, private communication; S. Dodelson, 
private communication), flnding that agreement between the 
models is better than ~ 3% (1.5% in temperature) over the 
range £ = 2 to 2000. 

We flt to the data using the 'eigenmodes of the signal-to- 
noise' approach discussed in Bunn & Sugiyama (1994), and 
similar to the approach of Bond (1994a, 1994b). We expand 
the map in terms of a set of basis functions, and use the coeffi- 
cients of this expansion to compute likelihoods for each model 
we vyish to test. The basis functions are chosen to have the 
maximum possible rejection power for incorrect models. (Note 
that we do not compute new eigenmodes for each new model 
we test: the same functions are used throughout.) We keep 
only the 400 most signiflcant modes in the expansion; retain- 
ing more modes does not signiflcantly increase the rejection 
power (Bunn & Sugiyama 1994). We use the weighted aver- 
age of the 31, 53 & 90 GHz two-year COBE maps, with a 
galactic cut at \b\ < 20°, and with the monopole and dipole 
removed. We have checked that other reasonable map combi- 
nations [e.g., using only the 53 and 90 GHz maps) give consis- 
tent answers. Our analysis is insensitive to the monopole and 
dipole in the maps, since we marginalize over these quantities; 
nonetheless, we remove a best-flt monopole and dipole from 
the maps. We have also performed Monte Carlo simulations 
to test whether this technique gives an unbiased estimate of 
{Q). We generated simulated data sets for an HZ power spec- 
trum with {Q) = 21/iK, and for an sCDM power spectrum 
with {Q) = 20/iK. We recovered a mean {Q) which differed by 
less than 0.1 /xK. 

We have applied this technique to determine the likelihood 
function for {Q) for a number of models. For sCDM we flnd 
{Q) = 19.9 ± 1.5/iK with the quadrupole included, and {Q) = 
20.3 ± 1.5 /xK with the quadrupole removed. In terms of power 
(Q2^ =403±61(/iK)2, with the quadrupole and (Q^) =418 ± 
64(/iK)2 without. For a tilted CDM model (with no tensors), 
the normalization is 

(Q)=21exp[0.69(l-re)]/iK (1) 

to an accuracy better than 0.5% for 0.5 < re < 2, while best flt is 
((Q),re) = (17.1/iK,1.3) [and ((<?), re) = (18.6/iK, 1.2) without 
the quadrupole]. The value of Q also depends weakly on Qb 
and is virtually independent of h, as indicated in the Table. 
Both Qb and h additionally affect the transfer function and 
hence the smaller scale porrnalizations. Our results are about 
1/iK higher than those in Gorski et al. (1994a). This appears 
to be due to the fact that we use the publically available sky 
maps rather than those used by the COBE group: when the 
Gorski et al. analysis is performed on the publically available 
maps, the results agree with ours to with 0.5 /xK. (Gorski, 
private communication; Gorski et al. 1994b). In any case, the 
discrepancy is well within the uncertainties. 

Including tensor modes reduces the goodness of flt to the 
data, since it increases the predicted quadrupole. With 
C2 /C2 = 7(1 — re) we flnd that the best-fltting values of 
(Q) = {{Q)l + (Q)^)i/2 are 22.9, 25.2 and 27.5/iK for CDM 
models with re = 0.9, 0.8 and 0.7, respectively. 

The 10° VARIANCE 

Given any set of C^'s one can predict the rms fluctuation, 
smoothed over 10°, through 

1 °° 

{a^{10°)) = -^{2l+l)C,W,, (2) 

t=2 
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where is the window function (Wright et al. 1994a). Thus 
it is possible to use it(10°) as a means of quoting the normal- 
ization. Using a simple 10° FWHM Gaussian for the window 
function leads to it(10°) = 1.99 (Q) for a pure re = 1 model. 
Using the window function and correction for the smoothing 
of Wright et al. (1994a) one has cr/{Q) = 1.91; the correspond- 
ing ratio for sCDM is 1.95. The data give us cr/{Q) = 1.63, 
which we find in less than 7% of a set of simulated COBE skies 
made with the sCDM power spectrum. Using re = 1.5 CDM 
instead of sCDM alleviates the problem: the predicted cr/{Q) 
is 2.32, while the data give a/{Q) = 2.31. 

However, if we remove the quadrupole from the data, then 
it(10°) and {Q) are quite consistent with each other. The ratio 
found in the data after removing the quadrupole is cr/{Q) = 
1.54, which agrees well with the prediction of 1.69 for sCDM. 
The lesson here is that it is better to normalize a model directly 
to the data, if it is available, rather than using the quoted 
10° variance. If you do use it(10°) then you certainly need to 
include the eff'ects of the beam of Wright et al. (1994a), which 
lowers cr/{Q) by ~ 5%. The remaining 'discrepancy' is due to 
the quadrupole being significantly lower than expected for fiat 
models (see Bennett et al. 1994). 

Other ways to quote the normalization. 

Let us take the spectrum of primordial fiuctuations to be a 
power law in comoving wavenumber k. If we assume that temp- 
erature fiuctuations at large scales arise from the SW effect, 
AT/T = 1/3$, the spectrum of radiation fiuctuations looks like 



P{k) = A{krioy 



(3) 



Here rjo ~ 3to = 2Hq^ (for f2o = 1) is the conformal time 
today, with scale factor normalized to unity and A is one way 
of quoting the amplitude for scalar perturbations. It is related 
to Sff, the dimensionless amplitude of matter fiuctuations at 
horizon crossing, through s^^- = (4/7r)A. 

Assuming the SW result above, we can write the average 
over universes of the moments of the temperature anisotropy 



C, = l^'-K^A 



r(3-re)r(£+^) 

r2(4^)r(£+s^) 



(4) 



(see e.g. Peebles 1982, Abbott & Wise 1984, Bond & Efs- 
tathiou 1987, White et al. 1994). For the special case of re = 1, 
we have C^j A = 47r/3 and oc £(£+ 1). Often one quotes 
not C2 but 



1/2 



1/2 



(5) 



which is the convention used in, e.g., Smoot et al. (1992), and 

where we use To = 2.726K (Mather et al. 1994). However, as 

has been emphasized before (Bond 1993), the C^'s for sCDM 

can depart significantly from equation (4) even for £ ~ 10. In 

fact they are not well fit by the SW formula for any value of 

re (Fig. 1). Note that for re > 1 the low-£ multipoles are quite 

suppressed relative to £ ~ 10, which goes some way toward 

explaining why larger re are preferred by the COBE data. 

Another common normalization convention is to define the 
matter power spectrum as 



P{k) = Bk"'T''{k), 



(6) 



where the transfer function T[k) ~ 1 on large scales. This 
means that the dimensions of B will depend on re (and will 



be Length for re = 1). If the fiuctuations arise purely from 
the SW effect, then B and A are simply related, so for re = 1 
(White et al. 1994) 

P{k) = 2n^4AkTl{k) (7) 
-2.5 X 10i*'A(A;//iMpc-^)r2(A;) {h'^M^cf . (8) 

So for an re = 1 primordial spectrum B = I-k'^tj^A = 
(67r^/5)?/Q {Q'^)/Tq. The assumption of pure SW fiuctuations 
is not exact (see Fig. 1), although when we normalize to the 
quadrupole it is approximately correct. A Boltzmann calcula- 
tion in which the photon and dark matter perturbations are 
explicitly evolved in time shows this relation holds to an accu- 
racy of < 4%. There is better agreement for higher values of 
h, which moves last scattering further into the matter domi- 
nated regime. For example, if /i = 1 the SW formula (for C2) 
is good to < 1%. We have used the ratio of matter to ra- 
diation normalizations from the Boltzmann calculation where 
appropriate. 

For CDM one conventionally uses a parameterized trans- 



fer function, T{k) = [l + {ak+ {bkf/'^ + {ckfY] 



(Efs- 



tathiou 1990) with a = (6.4/r)/i-iMpc, h = (3/r)/i-iMpc, 
c = (1.7/r)/i~^Mpc, v = 1.13 and the shape parameter F ~ 
Q,Qh. We use a calculation of the transfer function for our spe- 
cific model [h = 0.5, Q,b = 0.05) which fits the above form with 
r = 0.48. We find that variations in the transfer function with 
Q,B can lead to 10% changes in small scale power. 

Large-scale fiows also provide a measure of the power spec- 
trum: the variance of the velocity field in spheres of radius r, 
Vj^s(r) can be expressed as an integral over the power spec- 
trum (e.g., Peebles 1993). This tends to probe scales similar to 
the degree-scale CMB experiments. Bertschinger et al. (1990) 
estimated the 3D velocity dispersion of galaxies within spheres 
of radius 40/i~^Mpc and 60/i~^Mpc, after smoothing with a 
Gaussian filter on 12/i~^Mpc scales. So we also quote the 
normalization in terms of the quantities V40 and Veo, corre- 
sponding to the above procedure. 

On smaller scales, associated with clusters of galaxies, one 
conventionally quotes J3(10/i~^Mpc) (e.g. Peebles 1993), and 
the 'bias' hp or variance of the density field in spheres of 
8/i~^Mpc radius Ug, defined through 



Jo 



(9) 



where r = 8/1 ^Mpc. The variance of galaxies, possibly biased 
relative to the matter (5gai = h8p), is roughly unity on a scale 
of 8/i~^Mpc (Davis & Peebles 1983, Loveday et al. 1992). For 
sCDM the COBE best-fit gives Ug ~ 1.3, i.e. a relatively un- 
biased model. However this depends on the values of h 
etc. that are adopted, as indicated in the Table. 

Conclusions 

The quality of the data used to test theories of large-scale 
structure has improved markedly in recent years. There is even 
some hope that the CMB will soon allow us to make quanti- 
tative estimates of important cosmological parameters (Bond 
et al. 1994, Scott et al. 1994). In particular, the COBE data 
are of sufiiciently high quality that some care is required to 
determine the proper COBE normalization for any particular 
theory. When normalizing a theory such as sCDM, it is not 
correct to assume a pure SW HZ power spectrum, or even an 
re = 1.15 SW power spectrum. For sCDM and several of its 
popular variants, the normalizations in the Table are to be 
preferred. 
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There is an apparent discrepancy between the normaliza- 
tions in the Table and the measured temperature anisotropy on 
an angular scale of 10°: the value of the rms determined from 
the sky maps is approximately 20% lower than the ensemble- 
average values (it(10°)) in the Table. The agreement is im- 
proved by using the non-Gaussian COBE beam. Even so, 
Monte Carlo simulations show that this discrepancy is signif- 
icant at a confidence level of ~ 93%, but that removing the 
quadrupole from both the data and the simulations greatly al- 
leviates the problem. Since the value of it(10°) is so sensitive 
to the quadrupole, and since there is some reason to suspect 



that the quadrupole may be contaminated, cr(10°) is not a 
good quantity to use in normalizing models. In any case, an 
estimate of {Q) derived from a detailed fit to a model is to be 
preferred over it(10°), since the fit to {Q) makes use of all of 
the information in the data. 

We would like to thank K. Gorski, W. Hu and J. Silk for 
helpful discussions. The COBE data sets were developed by 
the NASA Goddard Space Flight Center under the guidance 
of the COBE Science Working Group and were provided by 
the NSSDC. This work was supported by grants from NASA, 
the NSF and the TNRLC. 
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3.26 


8.16 


1.36 
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nB= 1% 


3.32 


8.34 


1.39 
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1.45 


= 10% 


3.21 


8.03 


1.34 


19.8 


40.4 


349 


433 


450 


1.22 


h = 0.3 


3.23 


8.05 


1.35 


19.8 


40.6 


298 


357 


149 


0.67 


tCDM (S) 


4.19 


4.26 


1.57 


21.4 


41.0 


305 


376 


336 


1.03 


tCDM (S+T) 


2.85 


2.90 


1.81 


22.9 


43.6 


252 


310 


229 


0.84 


MDM 


3.26 


8.16 


1.36 


19.9 


40.6 


370 


463 


329 


0.97 



Table. — COBE normalizarions for standard Cold Dark Marror (sCDM) and several of irs variants. The firsr rliree nornralizarions use our besr fir for {Q'^) 
(with typical 'Itr* uncertainties ~ 15%), while the others are proportional to {Q) (error ~ 7.5%). The first two models were normalized using equation (4) for 
the angular power spectrum, neglecting everything but the Sachs- Wolfe effect. The third model is sCDM. The next three models show the effect of varying 
the baryon fraction (at fixed h — 0.5) and the Hubble constant (with Q^^h^ — 0.0125). tCDM is a tilted CDM model, with n — 0.9, given for the case with 
scalars only (S) and with a component of tensors added (S+T). MDM is a mixed dark matter model containing 30% hot and 69% cold dark matter, taken 
from Holtzman (1989). The power normalizations A and B are defined in equations (3) and (6) respectively. is the value of the angular power spectrum 
at i = 2. We have emphasized that the derived values of Q and <t are expectation values - this is strictly true of the other quantities also. Here by (t(10°) 
we mean equation (2) smoothed with a 10° FWHM Gaussian. T^o and T^q are the predicted rms velocities in spheres of radius 60/i~^Mpc and 40/i~^Mpc 
after smoothing with a Gaussian of width 12/i~-'-Mpc. The J3(10) and (Tg normalizations are defined in the text. 



